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ABSTRACT 



Aims. We investigate intensity oscillations observed simultaneously in the quiet chromosphere and in the corona, above an enhanced 
network area at the boundary of an equatorial coronal hole. 

Methods. A Fourier analysis is applied to a sequence of images observed in the 171 A and 1600 A passbands of TRACE. Four inter- 
esting features above the magnetic network are further investigated by using a wavelet analysis. 

Results. Our results reveal that, in both the 171 A and 1600 A passbands, oscillations above the magnetic network show a lack of 
power at high frequencies (5.0-8.3 mHz), and a significant power at low (1.3-2.0 mHz) and intermediate frequencies (2.6-4.0 mHz). 
The global 5-min oscillation is clearly present in the 4 analyzed features when seen in the 1600 A passband, and is also found with 
enhanced power in feature 1 (leg of a large coronal loop) and feature 2 (legs of a coronal bright point loop) when seen in the 171 A 
passband. Two features above an enhanced network element (feature 3 and feature 4) show repeated propagating behaviors with a 
dominant period of 10 minutes and 5 minutes, respectively. 

Conclusions. We suggest these oscillations are likely to be slow magneto-acoustic waves propagating along inclined magnetic 
field lines, from the lower solar atmosphere into the corona. The energy flux carried by these waves is estimated of the order of 
40 erg cnT 2 s for the 171 A passband and is far lower than the energy required to heat the quiet corona. For the 1600 A passband, 



the energy flux is about 1.4 x 10 erg cm 



. which is about one third of the required energy budget for the chromosphere. 
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1. Introduction 

The solar atmosphere is highly structured by the magnetic 
field. The chromospheric network, which is the upward exten- 
sion of the supergranular boundaries above the photosphere, 
is characterize d by clusters of magnetic flux concentrations 
(Gabriel 1976 ). Part of t he network flux e xpands with height 
(Patsoura kos et al.lll999t ITian et al.ll20 08a) and opens into the 
corona with a funnel shape, while the rest of the network con- 
sists of a dens e population of low-lying loo ps with lengths less 
than 1 4 km (iDowdy etalJll986t iPeteifcoOll) . 

The observation of oscillations in the chromosphere and 
corona can help us understand the magnetic structure of the so- 
lar atmosphere, and provides valuable insight into the unresolved 
coronal heating problem. 

It has been shown that oscillations with different frequencies 
(0.6 mHz-2.7 mHz) are p resent in the chro mospheric network 
Jcurdt and Heinzell |1998[ iLites et all 1 19931: ICurdt and Heinzell 
Il998t ICauzzie t al. 2000). McAteer et all (120021) found waves 
with a period of 4-15 minutes (1-4 mHz) in the central portions 
of network bright points, and suggest that these waves are possi- 
bly magneto-acoustic or magneto-gravity modes. More recently, 
McA teer et al.l (120041) found that the most frequent network os- 
cillation has a period of 283 s, with a lifetime of 2-3 cycles in 
four TRACE (Transition Region and Coronal Explorer) pass- 
bands centered at 1700, 1600, 1216, and 1550 A. The power 
of the intensity fluctuations in the 1600 A passband (UV con- 
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tinuum at 4 - 10 x 10 3 K plus Fe n at 1.3 X 10 4 K) which 
is formed at the temperature minimum (McAt eer et al]|2004l) . 
is found to peak at 3-5 mHz and considered as a mixture of 
photosphe ric 5-min oscillations and chromospheric 3-min os- 
cillations dFossum and Carlssonl2005l) . The so-called "magnetic 
shadows", which surround magnetic network elements and show 
a lack of oscillatory power in the period range of 2-3 min- 
utes, have received much attentio n ( Mcintosh and Ju dge 2001; 
Krijg er et al J 1200 it IVecchio et al.ll2007l) and have been consid- 
ered to play an important role in the trapping of high -frequency 
magneto-acoustic oscillations ( Sriv astava et al.ll2008l) . 

Oscillations with different periods have also been found in 
different regions above the chromosphere. In the active region, 
3-min oscillations have be en found in the transition region and 
corona in sunspot regions (lFludralll999t iBrvnildsen etd]fT999l 
iDe Moorte]| l2002). while 5-min o scillations are mo re likely to 
be found i n "non-sunspot" loops dDe Moorte]||2002l) and above 
AR plage (De Pontieu 2003, 2005). These oscillations are sug- 
gested to be due to upward-propagating slow magneto-acoustic 
waves. The 5-min oscillation is likely to be guided along in- 
clined magnetic flux tubes, which can decrease the acoustic cut- 
off frequency to allow the low-frequency p hotospheric oscil- 
lations propagate into the outer atmosphere dbe Pontieul r2004; 
iMcIntosh and Jefferiei2006tlHansteen et al]|2006l) . In the open- 
field regions, quasi-periodic perturbations with periods of 10- 
15 minutes dDeForest and Gurmanl 1 1998b and 10-25 minutes 
(Ban eriee et alj|20001) have been found in polar plumes. These 
low-frequency perturbations are explained as compressive waves 
(such as sound waves or slow-mode magneto-acoustic waves) 
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propagating along the plumes (Ofmanet al. 1999). In spicules, 
the 5-min oscillation is also present (|Xia et al.ll2005l) . and might 
be driven by p-modes (|Pe Pontieull2004l) . By measuring phase 
delays between intensity osci llations and between velocity oscil- 
lations of different line pairs. lO'Shea et al. 1 (120071) concluded that 
propagating slow magneto-acoustic waves are present in coro- 
nal holes and that they occur preferentially in bright regions that 
are associated with magnetic field concentrations in the form of 
loops or bright points. Periodic oscillations (7-64 minutes) in 
coronal bright points have also bee n fou nd bylUgarte-Urra et alj 
j2004al).IU garte-Urraet all d2004bl). andlTian et alj J2008bl) . By 
using the new XRT data, iKarivappa and Varghesd (120081) re- 
cently found that X-Ray bright point emission shows several 
significant peaks at different frequencies corresponding to time 
scales that range from a few minutes to hours. 

In this contribution, we investigate intensity oscillations ob- 
served simultaneously in the quiet chromosphere and in the 
corona, above an enhanced network area at the boundary of an 
equatorial coronal hole. 



2. Observations and data analysis 

The data set analyzed here contains a sequence of images ob- 
served from 07:40 to 08:20 UTC on November 08 in 1999 in 
two passbands (1600 A and 171 A) of the TRACE instrument. 
The emissions of the 1600 A and 171 A passbands are formed 
in the chromosphere (at temperature minimum) and in the lower 
corona, respectively. The sequence of 171 A only lasted from 
07:40 to 08:12. The pixel size is 0.5" for the 1600 A images, 
and 1.0" for the 171 A images. Standard software for calibrat- 
ing and correcting the TRACE data was applied to this data set, 
including removal of cosmic rays, subtraction of the dark cur- 
rent, normalization of the counts, and so on. We used the cross 
correlation technique to do the coalignment of the images in the 
sequence. 

We extracted a sequence of sub-images with a size of 130" x 
95" from the original data set. The sub-images were taken 
around the boundary of an equatorial coronal hole above an en- 
hanced network area. The average intensity images of the stud- 
ied region in the 1600 A and 171 A passbands can be found in 
Fig.|2j\ and E, respectively. The chromospheric network pattern 
is clearly seen in the 1600 A passband. The coronal hole and the 
surrounding quiet Sun can be easily discerned in the image of the 
171 A passband. Fig.Q]shows the magnetogram corresponding 
to the studied region. 

The cadence of our TRACE data was about 30 s. At a few 
moments there were irregular gaps. This problem was over- 
come by applying a linear interpolation with the help of the 
INTERPOL function in the IDL software. We then created a se- 
quence of running difference images, by subtracting each image 
from the image taken 30 s earlier. This sequence of running dif- 
ference images was further analyzed to search for periodic sig- 
natures. 

We applied a standard Fourier analysis to the sequence of 
running difference images of each passband. At each spatial 
pixel, we obtained a frequency-power curve. Then we added 
the total power respectively in three different frequency ranges, 
and calculated the percentage of this value relative to the total 
power in the full frequency range. In this way we obtained a rel- 
ative power map in each frequency range. The three frequency 
ranges are 1.3-2.0 mHz (500-780 s, low frequencies), 2.6-4.0 
mHz (250-384 s, intermediate frequencies), and 5.0-8.3 mHz 
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Fig. 1. MDI magnetogram of the studied region taken at 
08:03:30 UTC on November 08 in 1999. The white contours cor- 
respond to bright lanes seen in the averaged image of the 1600 A 
passband. 



(120-200 s, high frequencies). The power maps can be found 
in Fig. [2] 

The intensity images and maps of Fourier power reveal sev- 
eral interesting features. We chose four features for a more de- 
tailed analysis. They are shown in Fig. [2] and outlined in white. 
By checking the original TRACE images and the full-disk MDI 
magnetogram, we found that feature 1 might be a leg of a large 
coronal loop, feature 2 is a coronal bright point. Feature 3 and 
feature 4 are above an enhanced network element and their en- 
hanced power forms an elongated shape at low and intermediate 
frequencies, respectively. 

We averaged the intensity in the outlined rectangular area for 
each feature, and created the corresponding running difference 
light curve. These curves are shown in the upper panels of Fig. |4] 
to Fig. |7] Then we applied a Fourier analysis and a wavelet anal- 
ysis to each running difference light curve. The Fourier spectra 
are shown in Fig. [3] The solid and dashed curves correspond to 
the power of the 1600 A and 171 A passband, respectively. 

The localized nature of the wavelet transform allows us to 
study the duration of any statistically sign ificant oscillations as 
well as their period dBaneriee et alJuOOOl) . By decomposing a 
time series into time-frequency space, one is able to determine 
both the domi nant modes of variability a nd how those modes 
vary in time dTorrence and Compol 1 1998b . Here we also per- 
formed a wavelet analysis for the running difference light curves 
of the four features, in order to find the most reliable periods. 
We chose the Morlet wavelet function, defined as a sine wave 
modulated by a Gaussian window, for our analysis. The wavelet 
transform suffers from edge effect at both ends of the time se- 
ries. And this effect is important in regions defined as the "cone 
of influence" (COI). To check whether the periodic signatures 
present in the wavelet spectrum are real or not, we have to per- 
form a significance test. Here we chose a confidence level of 
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Fig. 2. Panel A-D: The average intensity image of the 1600 A passband, and maps of Fourier power in frequency ranges of 1 .3-2.0 
mHz, 2.6-4.0 mHz, and 5.0-8.3 mHz. Panel E-H: similar to Panel A-D, but for the 171 A passband. 
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Fig. 3. The Fourier spectra for the four features. The solid and dashed curves correspond to the power of the 1600 A and 171 A 
passbands and their scales are given on the left-hand and right-hand y-axes, respectively. 



95%. The description of t he wavelet function, COI, an d signifi- 
cance test can be found in lTorrence an d Compo d 1998b . 

The wavelet power spectra of the four features are shown 
in Fig.|4]to Fig. [7] Cross-hatched regions indicate the "cone of 
influence". The darker parts represent higher power, and the con- 



tours correspond to the 95% confidence level. We also averaged 
each wavelet power spectrum over time and plotted the global 
wavelet spectrum. 
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(a) TRACE 1600 intensity (a) TRACE 171 intensity 
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Fig. 4. Wavelet power spectra of the 1600 A (left) and 17 1 A (right) passbands for feature 1 . (a) The running difference light curve, 
(b) Time/period variation of the wavelet power spectrum. Cross-hatched regions indicate the "cone of influence". The darker parts 
represent higher power, and the contours correspond to the 95% confidence level, (c) Global wavelet. 



(a) TRACE 1600 intensity (a) TRACE 171 intensity 
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Fig. 5. Same as Fig.|4]but for feature 2. 



3. Results and discussion 

The main motivation of this work is to search for the signature 
of oscillations at the boundary of an equatorial coronal hole, and 
to study oscillations observed almost simultaneously in the chro- 
mosphere and corona above the magnetic network. The Fourier 
and wavelet analyses show that obvious oscillations really exist 
in our data. 

Fig. [2}3-D show the maps of Fourier power in different fre- 
quency ranges in the 1600 A passband. It is very clear that the 
network oscillations reveal a lack of power at high frequen- 
cies, and a significant power at low and intermediate frequen- 
cies. The reason why the power at high frequencies in the inter- 
network region is stronger than that in the network lane might 
be a result of the network magnetic structure. It is known that 
part of the network flux expands with height and opens into the 
coronadDowdv et al.lfl 986t lPeteiil200 It: iPatsourakos et al.|[l999t 
iTian etalj i 2008a). Thus, photospheric oscillations with frequen- 
cies above the acoustic cutoff (about 5.5 mHz) can propagate 



into the chromosphere, where the network fluxes begin to stretch 
into the region above t he internetwork in the shape of canopies 
(ISrivastava et al.ll2008l) . Below the canopy field, magnetic shad- 
ows with reduced oscillatory power at high frequencies are 
present within and immediately around the magnetic network 
(Vecch io et al 1 12001) . But the significant power (seen in Ca n 
854.2 nm) directly above the magnetic network in lVecchio et alj 
(2007]) is not very clear in our data set, which might be due to 
the different formation height between Ca n and the 1600 A pass- 
band of TRACE, or due to the relatively low spatial resolution 
of our data set. 



The result of a significant power at low and intermediate fre- 
quencies might be explained as a leakage of low-frequency pho- 
tospheric oscillations into the chromosphere through "magneto- 
acoustic portals" which are positioned within the network. It 
might pro vide a significant energy source to heat th e quiet chro- 
mosphere (iJefferies et al.ll2006t IVecchio et al.ll2007l) . 



H. Tian and L.-D. Xia: Network oscillations at the boundary of an equatorial coronal hole 5 



(a) TRACE 1600 intensity (a) TRACE 171 intensity 
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Fig. 6. Same as Fig.[4]but for feature 3. 




As for the power maps of the 171 A passband, we can make 
a similar conclusion to that of the 1600 A passband. But the 
significant powers at low and intermediate frequencies are more 
constricted and isolated in the network lane. Also, the regions 
with reduced power at high frequencies are not so prominent 
and are more diffused. It might be a natural result of the diffusion 
and expansion of the magnetic structure from the chromosphere 
to the corona. 

As we stated above, feature 1 corresponds to a leg of a large 
coronal loop system. Fig. [2jT and Fig. |2]j reveal a significant 
power at intermediate frequencies in this loop leg, both in the 
1600 A and 171 A passbands. The period corresponding to these 
frequencies is around 5 minutes. Although o scillations with peri - 
ods 1 80-420 s have already been reported bv lDe Moortel (|2000), 
it is still very interesting to note that in feature 1 the two bright 
dots seen in Fig. [2{j correspond exactly to the two parts of the 
loop leg seen in Fig.[2ji. This result indicates that in a loop sys- 
tem, the strongest power of the 5-mi n oscillation locate s at the 
lower part of the leg. As discussed in lDe Moortell d2000l) . the 5- 
min oscillation might be due to the propagating slow magneto- 



acoustic waves. From Fig. [3] and Fig. |4] we find that the power 
peak around 5 minutes is clearly present at most of the time. 
There is a weak peak at 10 minutes for the Fourier spectra in 
both the 1600 A and 171 A passbands, but not significant in the 
wavelet power spectra. 

Feature 2 corresponds to a coronal bright point, which 
is characterized by an enhanc ed coronal emissi o n assoc iated 
with bipolar magnetic fi eld. lUgarte-Urra et al.l (2004a) and 
Ugar te-Urra et al.l d2004bl) studied periodic oscillations of coro- 
nal bright points, and found oscillations with a period of 400- 
1100 s. L onger period (8-64 minutes) oscillations were also de- 
tected by Tia n et all (l2008bl) . Here from Fig. [3] and Fig. [5] we 
find that the oscillatory power of this bright point has a major 
peak around 5 minutes, and a second peak around 3 minutes. If 
we compare Fig. [2ji and Fig. |2]j, we can find that the strongest 
power of feature 2 corresponds to the two parts of the bright 
emission. These two parts might be the two legs of the loop sys- 
tem associated with the bright point. 

Feature 3 locates above the outer part of an enhanced net- 
work element. The most interesting point is that the enhanced 
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Fourier power at low frequencies has an obvious elongated shape 
in the 171 A passband. The movie of the image sequence clearly 
reveals that a periodic motion is present in this region, and the di- 
rection of this motion is exactly along this elongated shape. The 
Fourier and wavelet analyses both reveal a 10-min oscillation in 
the 17 1 A passband. This oscillation is present in the entire dura- 
tion of the time series, which can be seen in Fig. [6] While in the 
1600 A passband, the Fourier spectrum shows three peaks at 3, 5, 
and 10 minutes. However, only the 3-min and 5-min oscillations 
are considered as significant in the wavelet power spectrum. It 
has been known that strongly inclined magnetic fie lds can sig- 
nificantly decrease the acoustic cu t -off frequency ( De Po ntieu 
l2004t iMcIntosh and Jefferiesl 120061; lHansteen et alll2006l) . and 
thus allow high-frequency photospheric oscillations (above the 
acoustic cut-off frequency 5.5 mHz) to propagate into the chro- 
mosphere and corona. Fig. [2]2 gives us the impression that the 
magnetic field of feature 3 seems to be strongly inclined, and 
thus is very likely to guide the high-frequency oscillations from 
the lower part of the solar atmosphere to the corona. The weak 
peak at 10 minutes in the Fourier spectrum of the 1600 A pass- 
band, although not prominent in the wavelet spectrum, might 
still be the counterpart of the dominant 10-min oscillation seen 
in the 171 A passband. 

As for feature 4, which is above the same enhanced network 
element as feature 3 but in a different location, the power peak 
around 5 minutes is also clear in both the Fourier and the wavelet 
spectra for each bandpass. From Fig. [7] we can see that in the 
1600 A passband the period is initially around 4 minutes, then 
increases to about 6 minutes in the second half of the duration. In 
the 171 A passband, a period of 5-8 minutes is present through 
almost the entire duration. The strong power of the 5-min oscil- 
lation also seems to form an elongated shape in the 171 A pass- 
band shown in Fig.[2]j. A periodic motion is also present along 
this elongated direction when seen in the movie of the image 
sequence. 

The propagation speeds can be estimated if we plot the run- 
ning difference along the long side of the bars shown in the upper 
panels of Fig. [8] which indicate the directions of the propagat- 
ing 10-min oscillation in feature 3 and 5-min oscillation in fea- 
ture 4 in the 171 A passband. These plots are demonstrated in 
the low panels of Fig. [8] The slope of the dashed line provides 
an estimate of the propagating speed, which is about 32 km/s 
for the 10-min oscillation, and 58 km/s for the 5-min oscil- 
lation. These values are lower than the propa gating speed of 
the slow magneto-acoustic waves derived by De Moortel (2000) 
and|DeForest and Gurmanl (fl998). which are approximately 70- 
165 km/s at a bright loop-footpoint and 75-150 km/s in polar 
plumes, respectively. However, taking into account the line-of- 
sight effect of the measured propagating speed, our values are 
still of th e order of the cor onal sound speed, which is about 
150 km/s (De Moortel 2000). Since Alfven oscillations are es- 
sentially velocity perturbations and do not reveal themselves in 
intensity fluctuations, our periodic propagating signatures are 
most likely to be slow magneto-acoustic waves. 

Here we will roughly estimate the energy fl ux carried by 
these waves. We simply take the relationship in ISakurai et alj 
(120021) : 



Fourier Power(l .3 — 2.0 mHz) 



Fourier Power(2. 6 — 4.0 mHz) 



IT ? I 



si _ 2£_ = 2— 

/ p C s ' 



(1) 



where 6I/I, p' /p, 6v and C s are intensity fluctuation, density 
fluctuation, wave velocity amplitude, and sound speed, respec- 
tively. In our data, we found a value of about 0.015 for 6I/I. 
Given the sound speed of 150 km/s at the formation temperature 
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Fig. 8. Upper: the same as Fig.^ and Fig.|2]3. The bars indicate 
the propagating directions of the oscillations. Lower: plots of 
running difference along the long sides of the bars shown in the 
upper panels. The dashed lines are used to calculate the speed of 
propagating oscillations. 



of the 171 A passband, we can calculate the wave v elocity am- 
plitude as 1.05 km/s. Following Ofma net al.l d 19991) . the energy 
flux carried by slow magneto-acoustic waves can be calculated 
in the following way: 



F = Q.5p(6v) 2 C s , 



(2) 



where F and p represent ener gy flux and density , respectively. 
By using p = 5 x 10~ 16 g crrT 3 dDe Moortelf2 000) and substitut- 
ing the values of p and 6v, we obtained an energy flux of about 
40 erg cirT 2 s _1 for the propagating oscillations. These values 
are much lower than the required heat input to the quiet corona, 
which is about 3 x 10 5 erg crrT 2 s _I in With broe and Novesl 
(119771) . Thus, these waves are not appropriate candidates for 
coronal heating. 

For the chromospheric oscillation, if we take the value of 
0.04 for 51 1 1 and the sound speed of 16 km/s at the formation 
temperature of the 1600 A passband, we can calculate the wave 
velocity amplitude as 0.32 km/s by using Equation Q] Given a 
particl e number d ensity of 10 21 cm -3 at the temperature mini- 
mum (Hitei[2004), we can calculate the proton mass density p. 
By substituting the values of p, 5v, and C s in the chromosphere 
to Equation|2j we obtain a value of 1.368 x 10 6 erg cirT 2 s for 
the energy flux carried by the waves at the emission height of the 
1600 A passband. T his value is very close to that estimated by 
Jefferies et al. (2006) and is almost one third of the required en- 
ergy budget fo r the chromosphere. Thus, our results support the 
conclusion in iJefferies et al.l d2006l) that low-frequency magne- 
toacoustic waves provide a significant source of energy for bal- 
ancing the radiative losses of the ambient solar chromosphere. 

Finally, we have to point out that the relatively short length 
(low amount) of data restricts our discussion and lowers the 
solidness of the results. In order to fully understand network 
oscillations at boundaries of equatorial coronal holes, deep in- 
vestigation to more data sets is required in the future. 
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4. Summary and conclusion 

With the help of Fourier and wavelet analyses, we studied in- 
tensity oscillations observed simultaneously in the quiet chro- 
mosphere and corona, above an enhanced network area at the 
boundary of an equatorial coronal hole. 

Images of Fourier power reveal that, oscillations above the 
magnetic network show a lack of power at high frequencies (5.0- 
8.3 mHz), and a significant power at low (1.3-2.0 mHz) and 
intermediate frequencies (2.6-4.0 mHz) in both the 171 A and 
1600 A passbands. The former result suggests that "magnetic 
shadows" not only exist in the chromosphere, but also extend 
into the lower corona. The latter result supports the concept of 
"magneto-acoustic portals" within the network through which 
low-frequency photospheric oscillations can propagate into the 
chromosphere and corona. 

We also studied 4 interesting features in more detail. The 
global 5-min oscillation is clearly present in all of the 4 ana- 
lyzed features when seen in the 1600 A passband, and is also 
found with enhanced power in feature 1 (leg of a large coronal 
loop) and feature 2 (legs of a coronal bright point loop) when 
seen in the 171 A passband. Two features above an enhanced 
network element (feature 3 and feature 4) show repeated prop- 
agating behaviors with a dominant period of 10 minutes and 
5 minutes, respectively. The derived values of the propagating 
speed are of the order of the coronal sound speed. These veloc- 
ities, together with the compressional nature of the oscillation, 
suggest that our periodic propagating signatures are most likely 
to be slow magneto-acoustic waves. 

We calculated the energy flux carried by these waves and 
found a value of about 40 erg crrT 2 s _1 for the 171 A passband, 
which is only a small fraction of the total energy required to heat 
the quiet corona. However, for the 1600 A passband, the energy 
flux is about 1.368 x 10 6 erg crrT 2 s _1 , which is of the order of 
the required energy budget for the chromosphere. 
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